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Abstract — We study the performance of random linear network 
coding for time division duplexing channels with Poisson arrivals. 
We model the system as a bulk-service queue with variable bulk 
size. A full characterization for random linear network coding 
is provided for time division duplexing channels 11| by means of 
the moment generating function. We present numerical results 
for the mean number of packets in the queue and consider the 
effect of the range of allowable bulk sizes. We show that there 
exists an optimal choice of this range that minimizes the mean 
number of data packets in the queue. 

I. Introduction 

Reference IT] considered the use of network coding in 
channels in which time division duplexing is necessary, i.e. 
a node can only transmit or receive, but not both at the same 
time. This type of channel is usually called half-duplex, but the 
term time division duplexing (TDD) was used to emphasize 
that the channel is not used in any pre-determined fashion, but 
instead may vary the amount of time allocated to transmit and 
receive. 

In particular. Reference fl] studied the problem of trans- 
mitting M data packets through a link using random linear 
network coding with the objective of minimizing the expected 
time to complete transmission of the M data packets. Refer- 
ence |2| focused on the problem of energy consumption of 
this scheme showing that there exists, under the minimum 
energy criterion, an optimal number of coded data packets 
to be transmitted back-to-back before stopping to wait for an 
acknowledgment (ACK). 

The assumption in References fl] and fT\ was that the 
source had M data packets in its buffer before starting trans- 
mission. In a more realistic network setting, this buffer may 
sometimes empty or contain fewer than M packets awaiting 
transmission. Then, the source node must choose to either wait 
for additional packets to arrive, or take those packets in the 
buffer and start performing random linear coding. 

The problem of queueing for network coding systems has 
been considered previously to account for burstiness or losses. 
Some of this work considered the case in which feedback is 
available, e.g. |3|. References |4| and |5| studied a system 
with random linear coding, slotted time, and a Bernoulli arrival 
process. However, previous work has not considered timing or 
TDD constraints. Our work considers the problem in which 
the channel is TDD, where the time is not slotted and, more 
importantly, the service time depends on the size of the bulk. 



We study these scenarios for random linear network coding 
(RLNC) for TDD channels when the data packets arrive 
randomly at the source node according to a Poisson process. 
This problem can be modeled as a bulk service queue with 
a general service process. Bulk or batch service queues have 
been studied widely, e.g. |l6l. However, the service time for 
RLNC TDD will depend on the transmission time of both 
the coded packets and ACK packets, on the number of coded 
packets that are sent, and on the propagation time, as shown 
in |T|. This means that we have a bulk queue with a general 
service time, where the service time depends on the size of the 
bulk, which is not common in the existing studies. Reference 
|7| recently studied this problem with Poisson arrivals, calling 
it the M/G('"'^)/1 queue, where the size of the bulk can 
range between m and K. We build on this work to develop 
the queueing model of RLNC TDD. 

We provide numerical results for the mean queue size of 
the system for different choices of the arrival rate, and the 
pair (m, K). We show that the choice of the pair (m, A') can 
greatly impact the mean queue size of the system, suggesting 
that there is an optimal choice that minimizes the mean queue 
size. We show that fixing m = 1 and choosing K according 
to the arrival rate A is optimal in this sense. Also, we show 
that having a fixed batch size, i.e. m = K, is not the optimal 
configuration in general. 

The paper is organized as follows. In Section 2, we describe 
the system model and provide an expression for the moment 
generating function of the service time. In Section 3, we 
present the queueing model. In Section 4, numerical examples 
are provided. Conclusions are summarized in Section 5. 

II. System Model 

We consider each data packet to be of fixed-length, arriving 
to a source node through a Poisson process with rate A 
packets/s. Upon arrival, the data packet is placed in a buffer to 
await encoding and transmission to the receiver, as in Figure [T] 
The buffer forms a first-in-first-out (FIFO) queue. The service 
time of the queue is given by the time it takes to transmit 
a group of M packets taken from the queue using random 
linear network coding for TDD channels 0]. The size of the 
group of packets is variable, where m < M < K. The pair 
(m, K) constitutes the range of the bulk size or number of 
packets taken to perform random linear network coding |8|. 
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Fig. 1. Queue model studied in this work. 
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Fig. 2. 11] Structure of coded data packet: a header of size h bits, n data 
bits, M coding coefficients of size g bits each. 
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If the buffer has fewer than m data packets, the system will 
wait until m packets arrive before providing service. If the 
buffer contains more than K packets, the system will service 
exactly K packets. Finally, if the buffer has M packets with 
m < M < K, then the system will service M packets. 
Note that the service time depends on the number of data 
packets taken from the queue at any time, i.e. the service time 
distribution is general but it depends on the size of the batch 
being transmitted. Thus, we can use the bulk queueing model 
M/G^™'-^)/! developed in |,7J to study the problem. 

We consider that once the sender in a link has M data 
packets, it wants to transmit them reliably at a given link 
data rate R [bps]. The channel is modeled as a packet erasure 
channel. Nodes can only transmit or receive, but not both at the 
same time. The sender uses random linear network coding to 
generate coded data packets, which means that each coded data 
packet contains a linear combination of the M data packets 
of n bits each, as well as the random coding coefficients used 
in the linear combination. Each coefficient is represented by g 
bits. For encoding over a field size q, we have that g — log2 q 
bits. A coded packet is preceded by an information header 
of size h bits. Thus, the total number of bits per packet is 
h + n + gM. Figure |2] shows the structure of each coded 
packet. 

As in [1|, the sender can transmit coded packets back-to- 
back before stopping to wait for the ACK packet. The ACK 
packet feeds back the number of degrees of freedom (dof), that 
are still required to decode successfully the M data packets. 
Transmission begins with the M information packets taken 
from the queue, which are encoded into Nm > M random 
linear coded packets, and transmitted. If all M packets are 
decoded successfully, the process is completed and a new 
batch of data packets can be serviced. Otherwise, the ACK 
informs the transmitter how many dofs are missing, say i. 
The transmitter then sends Ni coded packets, and so on, until 
all M packets have been decoded successfully. 

Figure |3] illustrates the time window allocated to the system 
to transmit iV^ coded packets. Each coded packet CP(l,i), 
CP{2, i), etc. is of duration Tp. The waiting time is 
chosen so as to accommodate the propagation delay and time 
to receive and ACK. The problem of choosing the optimal 
number of coded packets to be transmitted back-to-back 
when i dofs are required at the receiver in order to decode the 
information was studied in [1]. 

This process can be modelled as a Markov chain. A state 
is defined as the number of dofs required at the receiver to 
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Fig. 3. Network coding TDD scheme. 



decode successfully the M packets. Thus, the states range from 
M to 0. This is a Markov chain with M transient states and 
one recurrent state (state 0). Note that the time and energy 
spent in each state depends on the state itself, because ^ 
Nj,yi ^ j in general. The transition probabilities from state 
i to state j iPi^j) are defined in \ \\. 

A. Moment Generating Function 

Let us define the moment generating function of the com- 
pletion time when the Markov Chain starts at state n as 



t 



exp(st)Pr(T = t) 



(1) 



where Pt{T — t) is the probability of the completion time 
being t. Note that Mx^^s) is the moment generating function 
of the completion time when n data packets are taken by the 
source to be transmitted reliably to the receiver 

Using the Markov Chain structure of the problem, it can be 
shown that My„(s) can be re-stated as 



exp . u / 

\ i=l 



mn>lm„„x>0 mx>0 

(2) 

where T* is the deterministic time required to send iVj coded 
packets and wait for an ACK when the Markov chain is in state 
i, i.e. T* = NjTp + T^. The constant C„ captures the effect 
of returning to the same state repeatedly, while An captures 
the different paths that can be traversed without repetition of 
a state. 

The expression for Cn is 



The coefficient for An can be shown to obey a recursive 
expression of the form 
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1 when s e 5 and zero otherwise. 

Substituting expression ^ into (|2| we obtain the following 
recursive equation for the moment generating function 
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with MTflis) = 1. 

Finally, note that the same structure is valid for computing 
the energy needed to complete transmission. To do so, one 
would substitute by E\ and Mj-nis) by Af£; „(s). 

III. QuEUEiNG Model 

The system model discussed in the previous section corre- 
sponds to the bulk queueing model M/G*^™'^-'/! studied in 
Q. This bulk queueing model considers Poisson arrivals and 
a general service time that depends on the bulk size. 

The transition probability of the number of packets in the 
queue is given by |7| 



{seS} 1^ which provides an expression for n(z) in terms of its first 
K + 1 coefficients ttq , . . . , ttj^ ■ Determining these K + 1 
coefficients provides a full characterization of the stationary 
probabilities Q. Reference Q proves that ^(^)(z) — has 
n-l exactly K zeros satisfying |z| < 1 assuming that A^-^^z) has 

E^__^^Mj' i{s) (3) ^ radius of convergence greater than one. Denoting the roots 
j_Q as 1, zi, and assuming that they are different, note 

that the numerator of (|7| has to vanish for zi, z^-i which 
gives \xs K — 1 linear equations 
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where is the probability of k arrivals during a service of 
type j. 

Let us define 



k=0 



(4) 



We can use a similar analysis as that of Reference |7| to 
prove that 

A^^\z) = Mt,,{X{z ~ 1)) (5) 

and that 
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The system is stable if and only if A < Kfij^, where I/ij^k 
is the mean service time when the bulk size is M = K, where 

Let us denote by Il{z) — X^i^o"^*^* corresponding 
generating function of the stationary probabilities. Reference 
Q showed that n(z) can be expressed as 
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for k = 1, ...K — 1. Also, the numerator vanishes trivially for 
z = 1 for both the numerator and the denominator in (|7|. 
We thus need one more linear equation. To obtain this we 
use I'Hospital's rule to exploit the fact that 11(1) = 1. This 
translates to 
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where we have used the fact that — » ^ ^ 

X/IJ'K K in general, the denominator will not become zero 
as z — > 1. In fact, if the system is stable this condition will 
be satisfied. 

The final linear equation to fully characterize n(z) given in 
Reference Q is 
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A. Queue of Finite Capacity 

The general solution requires the calculation of the roots of 
AW(z) — z , which can result in numerical inaccuracies in 
practice because has exponential terms. Also, calculating 
the roots is increasingly difficult when the decision variable 
K assumes a larger value. For these reasons, we simplify 
the problem considering that the system has a capacity of B 
packets waiting to be serviced, i.e. without considering those 
that are being transmitted. The transition probability for this 
case is 
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(7) where R{k, I) ^ 1 



Ek 
j=0 



AO 



0.]HL 



if. 

r' 'tf 
, I i,\ 
'i 1,1 

r M+ 
9, \ 




m = 5,3,1 



* » i n> » > » I II » 1 I t 1 » » » a » » » » 1 

10 15 20 25 30 



Queue Size (i) 

Fig. 4. Stationary distribution 7ri,j = 0,1,. for A = 30 packets/s, 
K = 5, B = 30, and different values of m. 



In order to compute the stationary distribution, it suffices 
to solve 7f = Ptt, with tt = [ttq, tti, tt^]"^, under the 
constraint that J2f=o'^i ^ ^■ 
B. Performance Analysis 

We will consider two metrics in order to study performance 
of the system. First, the mean queue size defined as 
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for the case in which the queue capacity is B. If there is no 
constraint on the capacity, we simply let i? ^ cx). 

The second metric is the mean batch size in steady state, 
which can take values {m,m + l, K}. Defining Z^,.^ as 
the batch size for a choice of (m, K), then 

m K-1 B 

E[Z(^rn,K)]='m'^'^i+ E in + K ^ TTi, (11) 
i=0 i=m+l i=K 



where X^i^iC = 1 ^ SiLo^ ^i' ^'^'^ there is no constraint 
on capacity, again we let _B — > oo. 

IV. Numerical Results 

This section provides numerical examples that show the 
performance of our network coding scheme for different 
settings of (m, K) and arrival rate A. We use the mean queue 
size defined in the previous section as our metric of interest. 
We use a high latency channel with packet erasure probability 
Pe — 0.2, propagation time of 12.5 ms, data packets of 
10,000 bits, g = 20 bits, a rate R = 1.5 Mbps, a header 
of /i = 80 bits, and the ACK packet has 100 bits. We assume 
that the number of coded packets to be sent back-to-back (iVj) 
are chosen to minimize the mean transmission time as in 

Figure |4] shows the stationary distribution for A = 30 pack- 
ets/s, K — 5, B — 30 and different values of m. This figure 
shows that the high probability states correspond to small 
number of packets in the queue. The probability of large queue 



sizes when the system operates in steady state is very low. 
Figure |4] also shows that for low values of to, the stationary 
distribution is concentrated in the lower values of the queue 
size i. As to increases, the stationary distribution spreads over 
larger values of the queue size i. 

Table |l] shows the mean queue size when A 1 packet/s 
under different configurations of the pair (to, if), with to < 
K. We observe that the mean queue size shows greater 
dependence on the value of to than on the value of K. For 
example, increasing the value of K when m = 1 shows 
little variation in the mean queue size, while increasing to 
with any value of K > 1 increases the mean queue size. In 
terms of minimizing the mean queue size for low values of 
A, this means that we should allow transmission of bulks of 
size 1 . This is the case because waiting for additional packets 
before transmitting is costly, considering that the time between 
packet arrivals might be larger than the mean service time 
for a single data packet. For the case of bulks of size 1, 
our scheme transmits several copies of the packet back-to- 
back before stopping for an ACK, which is similar to the idea 
presented in |9|. 

Table |ll] shows the mean queue size when A = 30 packet/s 
under different configurations of the pair {m,K), with to < 
K. For this A, we do not consider the setting m — K = 1 
because it is not stable for the infinite capacity case {B oo) 
and will present very high packet drops when the capacity B 
is finite. Again, there is an advantage of allowing to = 1 for 
the studied cases, in terms of reducing the mean queue size. 

Tables U and |ll] also show the mean batch size for A = 
1 packet/s and A = 30 packet/s, respectively. The mean batch 
size is biased by the value of m and K. However, it gives us 
some intuition about the operation of the system. For example. 
Table |l] shows that when m is too large with respect to the 
arrival rate, the mean batch size is close to to. This means that 
on the average the system services the batches much faster than 
the time it takes the queue to fill with to new packets. Thus, 
the system will be idle for long periods of time just waiting for 
the queue to fill to the required m. Only with small probability 
the batches will contain more than to packets. Of course, if 
m ^ K the batch size will always be to as seen in the tables. 

Let us consider the case of a fixed batch size, i.e. m — K. 
Tables [l] [ll] and Figure |5] show that the optimal choices are 
TO = X = 1, TO = if = 2, and to = X = 3 for 
A = 1 packet/s, A = 10 packet/s, and A = 30 packet/s, 
respectively. However, we notice that the mean queue size is 
larger than other configurations of (to, K) without the fixed 
batch size restriction. For example, the optimal fixed batch 
size configuration for A = 30 packet/s is 53 % larger than the 
{m,K) = (1,5) configuration presented in the table. Thus, 
we observe that having a fixed batch size is not the optimal 
configuration in general. 

V. Conclusion 

This paper provides a queueing model for random linear 
network coding scheme for time division duplexing channels 
with Poisson arrivals. The analysis considers that the size of 
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TABLE II 

Mean Queue Size and Mean Batch size for different (m, K) 

CONFIGURATIONS. THE PARAMETERS USED ARE A = 30 PACKET/S 
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the batch that is sent using random linear coding can be in a 
range of values, say (to, K). At the time of completing service 
to a batch, if the queue size is below the minimum allowable 
value of the batch to, the system will wait until the queue 
size becomes m. If the queue size is in the range of batch 
sizes, all data packets are serviced at that time. Finally, if the 
queue size is greater than the maximum allowable batch size 
K, the first K packets of the queue are serviced. We have used 
this analysis to study the mean queue size of the system and 
to choose the (to, K) pair that minimizes it under different 
A. The analysis is useful if we are interested in choosing the 
optimal (to, K) given a different objective function. 

We present a complete characterization of the problem of 
random linear network coding for time division duplexing 
presented in 11], by providing a recursive expression for the 
moment generating function of the service time. This moment 
generating function is also valid for the energy to complete 
transmission of a batch of data packets using the appropriate 
substitutions. 

Numerical results suggest that the mean queue size shows 
greater dependence on the value of minimum batch size to 
than on the value of the maximum batch size K for low values 
of A. In general, K determines the maximum serviceable A of 
the system, while m should be allowed to have small values. In 
our examples, to = 1 provided the best performance in terms 
of minimizing the mean queue size. Also, numerical results 
suggests that having a fixed batch size, i.e. to = K, is not the 
optimal configuration in general. 
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Fig. 5. Mean queue size for the fixed batch size case (m = K) with B = 30. 

Future research should consider choosing (to, K) to opti- 
mize the total transmission time and total transmission energy 
of a packet, i.e. from the time it enters the queue until it 
is serviced. Also, extensions of the principles proposed for 
one link to the general problem of wireless networks will be 
studied. These extensions are possible within the framework 
of random linear network coding. 
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